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Abstract 

Multi-phonon Raman scattering in semiconductor nanocrystals is treated 
taking into account both adiabatic and non-adiabatic phonon-assisted optical 
transitions. Because phonons of various symmetries are involved in scattering 
processes, there is a considerable enhancement of intensities of multi-phonon 
peaks in nanocrystal Raman spectra. Cases of strong and weak band mix- 
ing are considered in detail. In the first case, fundamental scattering takes 
place via internal electron-hole states and is participated by s- and d-phonons, 
while in the second case, when the intensity of the one-phonon Raman peak 
is strongly influenced by the interaction of an electron and of a hole with in- 
terface imperfections (e. g., with trapped charge), p-phonons are most active. 
Calculations of Raman scattering spectra for CdSe and PbS nanocrystals give 
a good quantitative agreement with recent experimental results. 
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I. INTRODUCTION 



The symmetry-breaking coupling of vibrational modes to degenerate electronic states 
described first by Jahn and Teller |l] has been manifested in a wide range of transport 
and optical phenomena in various electron-vibrational systems. The role of the Jahn- Teller 
effect in the optical transitions of semiconductor nanocrystals has been recently studied 
in Refs. fHH for the strong confinement regime, when the radius of a nanocrystal R is 
smaller than the Bohr radius of an exciton in the bulk R ex , and in Refs. MM for the weak 
confinement regime, when R > R ex . In the case of a strong confinement, when all the states 
of an exciton, including the ground state, are degenerate due to band mixing, the Jahn- 
Teller effect substantially influences the phonon-assisted optical transitions already in the 
framework of the dipole approximation As a result, the optical spectra of nanocrystals 
with R < R ex can be different from the Franck-Condon progressions typical of the adiabatic 
systems. 

Recently, multi-phonon peaks in Raman scattering spectra have been observed in CdSe, 
CdSe^Si-a, and PbS nanocrystals (see, for example, Refs. Previously, the theoretical 

analysis of multi-phonon Raman scattering in nanocrystals has been generally based on the 
adiabatic approximation due to the seminal works of Huang and Rhys || and Pekar []T0 



However, even taking into account band mixing, the values of the Huang-Rhys factor for the 
nanocrystals of radius R = 1 nm to R = 5 nm appear to be considerably smaller than unity. 
This result is at variance with those derived from the experiments cited above. In the present 
paper, the theoretical treatment is performed involving both adiabatic and non-adiabatic 
optical transitions. The interaction of phonons with an exciton in a degenerate state leads 
to internal non-adiabaticity of nanocrystals (the Jahn- Teller effect), whilst the existence of 
exciton levels separated by an energy comparable with that of the optical phonon results in 
external non-adiabaticity (the pseudo Jahn- Teller effect). The structural imperfections of the 
real interface between a nanocrystal and a host medium can also substantially influence the 
probabilities of phonon-assisted transitions, see e. g. Ref . [|11[] . Moreover, these imperfections 
(e. g., trapped charge, localized at the interface, as supposed in Refs. [jSUl2|) can increase the 
contribution of non-adiabatic transitions into the scattering probability, as shown below. 

In nanostructures, phonon spectra drastically differ from those in bulk materials 
PJT3|-p^o]| . Previously, several models of polar optical phonons and of the electron-phonon in- 



teraction have been developed for these structures. The dielectric continuum model which is 
proposed in Ref. JlTj and used for planar fi~8l-f22f , cylindrical and spherical |6j,^3|l structures, 
exploits only electrostatic boundary conditions. The alternative hydrodynamic model (see 
Refs. |p^-p7|) treats only mechanical boundary conditions. Distinctions between the two 
models consist in a choice of a complete orthogonal basis for the relative ionic displacement 
u(r). It has been shown in Refs. [^,[^J that physical quantities expressed as a sum over 
all phonon modes (scattering rates, polaron parameters etc.) slightly depend on a concrete 
basis. On the other hand, calculated Raman spectra are sensitive to mechanical boundary 
conditions. Various improvements of the dielectric continuum model [31], 31] are developed in 
order to obtain better agreement between macroscopic and microscopic p9| , |32]| description 
of lattice vibrations and with experimental data on Raman scattering [[33H 3511 ■ Within the 



continuum model of optical phonons in spatially confined systems (Refs. |T5| . p!6| , |36| -p^] ) , both 
electrostatic and mechanical boundary conditions are imposed on the ionic displacement. 
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The multimode dielectric model developed by us (see Refs. |39 , |40| ), in addition, explicitly 
takes into account a finite number of phonon degrees of freedom in a quantum dot. 

In Section 2, the set of phonon modes is determined using an effective dielectric function 
of a nanocrystal in a finite-dimensional basis for the field of the ionic displacements. In 
Section 3, multi-phonon resonant Raman scattering in spherical nanocrystals is considered 
using the obtained Hamiltonian of the electron-phonon interaction. The obtained analytical 
and numerical results are discussed and compared with experimental data in Section 4. 



II. OPTICAL PHONONS IN NANOCRYSTALS 

In the present investigation, longitudinal optical (LO) vibrations in spherical quantum 
dots are considered on the basis of the multimode dielectric model (see Refs. PP|^D| ). Polar 
optical phonons in a semiconductor nanocrystal of a radius R embedded into a polar medium 
can be described within the continuum approach (when R 3> a , where ao is the lattice 
constant) by the renormalized relative ionic displacement vector 

m M e >k (r) 

w fe ( r ) = j I 1 ) 

VofcWfc,TO\/£o [ £fe (°) _ £k (°°)] 

where k = 1 for the nanocrystal, k = 2 for the host medium. For the k-th medium, u/~ (r) is 
the relative ionic displacement, e* k is the effective charge of an ion, uOk,io is the frequency of 
the transverse optical phonon in the Brillouin zone center, Vofc is the volume of the elementary 
cell, Ek(oc) and £fe(0) are the high-frequency and static dielectric constants, respectively; e 
is the permittivity of vacuum. Dynamics of the non- dispersive polar optical phonons is 



determined by the Born-Huang equation |41] combined with Maxwell equations. In order 



to analyze the time evolution of w*. (r) taking into consideration the spatial dispersion, the 
generalized equation of motion is used: 

(^fc,TO - w fe ( r > u ) = ^£;,to\Ao [£k (0) - e k (oo)]E fc (r, uj) 

+ J T k (r - x) w k (x, uj) dx. (2) 

The right-hand side of Eq. (H) consists of two terms: (i) the long-range force expressed by 
means of the macroscopic electric field (r, uj) induced by the ionic displacement vector 
Wfc (r, uj); (ii) the dispersion term of the short-range force with the spatial dispersion tensor 
Tfc (r). The action range of this force can be estimated as a few lattice constants. 

The set of equations to determine the polar optical phonon dynamics in both the 
nanocrystal and the host medium consists of the equations (0), the electrostatic Maxwell 
equations, and the equation relating the electric field and the polarization to the electrostatic 
displacement for each medium. The form of the spatial dispersion tensor T^ (r) is chosen in 
such a way that the bulk LO (TO) phonon dispersion culo(to) (q) be reproduced when the 
equation of motion (|2|) is solved in bulk. Further on, anisotropy of the dispersion of bulk 
phonons is neglected: the bulk phonon frequencies depend only on the wave number q. 

For a nanocrystal, the aforesaid set of equations must be completed by electrostatic and 
mechanical boundary conditions. In the present work, we choose the following mechanical 
boundary conditions |p7j|39 
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w~j: = at the interfaces, k = 1,2. (3) 
They lead to a double hybridization of LO and interface modes. For planar multilayer 



structures, phonon spectra obtained in Ref. |37| when applying Eq. (^) are in excellent 



agreement with experimental Raman scattering data ||42| . 

Owing to the spherical symmetry of a quantum dot, the eigenvectors of phonon modes 
g/ m (r, u) correspond to a definite phonon angular momentum I and its z-projection m. The 
index J labels LO and TO phonon branches. We seek these eigenvectors as a superposition 
of orthogonal vector functions ff ms (r) which satisfy Eq. (^|): 

g/ m (r,^)=£f/iHf/ ms (r). (4) 

s 

This basis is subdivided into three subsets describing: 

(i) LO phonons with definite values of the angular momentum I, of its projection m, and of 
the radial quantum number s, 

^Ims ( r ) = 




-Ji 



R 



r(l) 
Im 
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(5) 



where n = r/r, ji (x) is a spherical Bessel function, b^ s is the s-th zero of its derivative j[ (x), 
Y^ 1 (n) (A = 0, ±1) is a spherical vector of a definite parity; 
(ii) TO phonons of the "electric" type, 



rTO,E 



Ims 



r) 



ai jS Jl + 1 [dig) 



■I' \ "IT") Y L 1} l") 



R 



+ 



R . b l>s r 



ai,,r 



Yffi (n) 



where is the s-th zero of j/ (x); 

(iii) TO phonons of the "magnetic" type, 



f TO,M 



Ims 



r) 



v/2 



(6) 



i? 3/2 Jw(a ilS ) J ' V -R 

Dimensionalities of the subsets are implicitly determined by the inequalities: 

bi >s < TrR/ao, for LO modes; 
ai >s < 7iR/ao, for TO modes, 



(7) 



(8) 



which express the fact that the "wavelength" of an optical phonon cannot be smaller than 
the double lattice constant. 
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Using the expansion (£|) for a joint solution of Eq. (0) and of the Maxwell equations with 
electrostatic boundary conditions, we arrive at the dispersion equation 



l + l 



0. 



(9) 



where e 2 (oS) is the dielectric function of the host medium. The function e\ (I, to) is deter- 
mined by the formula 



1 
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l,LO 



(0) 



l,TO 



(0)' 



where Q 



1,8 



e\{l,u) ex{oo) 
bi }S /R, and the coefficients xi,s are: 

21 



CO 



l,LO 



(Qi,s) 



00* 



(10) 



(ii) 



These coefficients possess the property Y^Li Xi,s — 1- 

We can interpret E\ (l,oo) as the effective dielectric function of the quantum dot, using a 
formal analogy of the dispersion equation @ with that of the dielectric continuum model (cf. 
Eq. (13) of Ref. ||). In the "non- dispersive" limit, when we set in Eq. ( |IUD u^lo (Ql,s) — 
^i,lo (0) , the function E\ (I, u>) becomes equal to the "bulk" dielectric function ei (to) , and 
Eq. (||) turns into the dispersion equation of the dielectric continuum model ||. The latter 
equation provides, in particular, interface phonon frequencies u)j t i(l) [j = 1,2; wi/ro (0) < 

LO jtl (I) < LOi >h0 (0)]. 

Basis vectors of phonon modes are denoted as gf (r), where the index rj labels the 
roots of Eq. (§) ui >f) at a definite value of the angular momentum I. Generally speaking, LO 
phonon modes cannot be subdivided into bulk-like and interface ones. Therefore, they can 
be considered as hybrids of both these types. The following explicit expression is obtained 
for the coefficients in Eq. (§): 



f/, LO 



Ls 



Ci {uji >v ) 



bl -1(1 + 1) col LO (Qi >s ) - col v 



(12) 



The transformation (| 
then found to be 



is unitary (see Ref. |39|). The normalization constant Ci (&>j,jj) is 



E 



[Ks -ni + i) 



LO 



l,LO 



(Qi,s 



— LO 



l,V 



-1/2 



(13) 



Finally, the phonon Hamiltonian takes the form: 

H L = J2 huJ 



while the Hamiltonian of the electron-phonon interaction is 

H e -L = £ [%a u + 7+0+) , 



(14) 



(15) 
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where the set of indexes v = (I, m, rf) labels the obtained LO modes. The amplitudes of the 
electron-phonon interaction can be explicitly expressed in terms of the unitary transforma- 
tion coefficients as follows: 



2Hu\, 



LO 



1/2 



R 



E 



\mu 1:LO 



V2 



1/2 



Y lrn (n) 



#.-1(1 + 1) 



ji {bi, s r/R) 

Jl (k,s) 



{l + l)e 2 (w t . 



lex + (l + l)s 2 (uj ltV ) \R 



(16) 



where a\ is the Frohlich electron-phonon coupling constant of a nanocrystal. The Hamilto- 
nian of the exciton-phonon interaction has the form similar to Eq. flT5| ) with the replacement 
of % by the exciton-phonon interaction amplitudes 



(3 V (r e , v h ) = % (r e ) - % (r h ) 



(17) 



It is worth noting that in CdSe and in PbS, where dispersion of bulk LO phonons is 
strong (see Refs. P^[43|fP]] ) , the LO phonon modes as derived from Eq. (|D|) appear to be 
hybrids of bulk-like and interface vibrations. 



III. RAMAN SPECTRA 



Within the context of the long-wavelength approximation, the interaction of an electron 
with an electromagnetic field is described by the operator V (t) = Vit~ %Vilt + Vge lUst , where 
the terms Vj and Vg correspond, respectively, to the absorption of a photon with the fre- 
quency Qj (incoming light) and to the emission of a photon with the frequency Qs (scattered 
light). The interaction amplitude Vj(S) is proportional to the projection of the electron dipole 
momentum operator d on the polarization vector e 1 ^ of the relevant wave: d 1 ^ = e 1 ^ ■ d. 
From the second-order perturbation theory, the transition probability between an initial \i) 
and a final |/) state is 



2tt 



E 



./ 



v 



m) (m 



V, 



(18) 



Here, u m i and Ufi are transition frequencies. For phonon-assisted Raman scattering in 
a semiconductor nanocrystal, \f), and \m) are quantum states of the exciton-phonon 
system. Both the initial and the final states contain no charge carriers (electrons or holes), 
so that these states are described by a direct product of a wave function of free phonons 
with that of the exciton vacuum. Intermediate states \m) are eigenstates of the Hamiltonian 



H = H ex + H ph + if, 



-phi 



(19) 



where H ex is the Hamiltonian of an exciton, H p h is the phonon Hamiltonian (|14|), and H ex _ p h 
is the Hamiltonian of the exciton-phonon interaction. 

For definite polarizations of the incoming and the scattered light, the scattering prob- 
ability is obtained by averaging Eq. fll8|) over the initial states and by summing over the 
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final ones. Since the phonon Hamiltonian Hi is quadratic and the Hamiltonian of the 
exciton-phonon interaction H ex -L is linear in phonon creation and annihilation operators, 
an analytical averaging over the equilibrium phonon ensemble can be performed in Eq. 
(Pf). As a result, we can express the shape of the Raman spectrum in terms of time-ordered 
operators averaged over exciton states only: 



oo oo 



W(fij,e/, fi s ,e s ) = J dt e ms ^ l)l j dr J da e -^+-)-^/(— ) 

-oo 



, t—lj„i a 



nlK^TsTjexp {$[$,$']} 



/'i 



/4 



(20) 



Here, lu^ are the frequencies and d^ s > = (/^tr^ |0) are the dipole matrix elements for a 
transition from the exciton vacuum |0) to the eigenstates of the Hamiltonian H ex . The 
"influence phase" of the phonon subsystem <3> f3, j3' is the following operator: 



h 2 



£ / ds / ds'T^ (f + a - X {s') Pi (s) 



ds / ds'T^ (s - s') X (s) ft (s') - ds ds'T* (s - s') P' v (s') # (s) . (21) 



In Eqs. ( P0|) and (^T|), primed (non-primed) exciton-phonon amplitudes are 
(anti) chronologically ordered operators, which act, respectively, on primed (non-primed) 
exciton states. The phonon Green's function 



T u (t) 



cosh [u (t - h/2k B T)} 
sinh (/ko>/2/c s T) 



(22) 



describes the phonon emission and absorption processes. 

Because electro n-phonon coupling is weak in nanocrystals, such as CdSe, CdSe^Si.^ and 
PbS, the /T-phonon scattering intensity, corresponding to a definite combinatorial frequency 



K 



w v ., can be analyzed to the leading (K-th) order in the electron-phonon coupling constant 



i=i 



The scattering intensity within this leading-term approach is then expressed through a 
squared modulus of the scattering amplitude: 



E 



d 1 (d s 



K 

n 

MO j=l Q 



Pv 



Hj-l 



(23) 



k=l 



where V ^ is the inverse lifetime of an exciton in the state \fj), while T u is the inverse lifetime 
of a phonon of the mode v. 

The following treatment is related to the regime of strong confinement R < R ex , where 
R ex is the exciton Bohr radius. In this case, the Coulomb electron-hole interaction can be 
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treated as a perturbation, and the concept of the exciton means the same as that of the 
electron- hole pair. 

If one considers states of the electron-hole pair in Eq. ( p3|) within the model of simple 
bands, then diagonal matrix elements of amplitudes of the electron-phonon interaction on 
the wave function of the ground electron-hole state vanish, what suppresses the one-phonon 
Raman scattering. This suppression can be removed by one of the following mechanisms: 
(i) the Coulomb electron- hole interaction ||46||, (ii) band mixing [J47]|, (iii) the influence of 



interface imperfections on the exciton wave functions pj| , (iv) non-adiabatic transitions 
with the participation of virtual phonons. The latter mechanism exists only beyond the 
leading-term approach. 

Two particular cases are of interest when one of these mechanisms dominates: (i) strong 
band mixing (CdSe in glass 0,01), (ii) weak band mixing (PbS in PVA ||) when one- 
phonon transitions are determined mainly by the scattering of an electron and of a hole by 
the potential due to the imperfections. It is worth recalling, that in the strong-confinement 
regime, the Coulomb interaction only slightly influences the electron-hole states and the 
Raman spectra. 

For a CdSe nanocrystal, exciton states are analyzed here using the spherical model of 
the exciton Hamiltonian K| supplemented by terms describing the electron-hole exchange 



interaction [49]: 



2m e 2m Q 9m v 

2 ^ 
+V C (r e , r h ) - -e exch a 3 5 (r e - r h ) (3 • j) , (24) 

where p e and f>h (r e and r^) are the electron and hole momenta (coordinate vectors); 71 
and 72 are the Luttinger parameters, m and m e are the bare electron mass and the electron 
band mass; P( 2 ) and are irreducible second-rank tensors of the momentum and of the 
spin-| angular momentum of a hole, respectively. The potential of the Coulomb attraction 
between an electron and a hole in a spherical nanocrystal Vc (r e , r^) is used from Ref. ||39|| . 
In Eq. (^), the last term, which is proportional to the scalar product of the electron (a) and 
hole (J) spin operators, describes the electron-hole exchange interaction characterized by the 
strength constant e exc h- In CdSe, this constant is equal to 320 meV The Coulomb and 
the exchange interactions between an electron and a hole are treated as perturbations when 
determining the exciton states. To the zeroth order in these perturbations, the exciton states 
\fi) are characterized by a definite electron spin projection a and a total angular momentum 
of a hole F with the ^-projection M: 

= a- K, M) = *ts < (M) , (25) 
where K = nS 3 / 2 , riPi/2, nP 3 /2, nP 5 / 2 , etc. The index n labels the solutions of the equations 



for the radial components of the hole wave function [[|7j . 

In contrast to CdSe, the conduction and valence bands in PbS [|5(J are both non- 
degenerate. As a result, band mixing in PbS nanocrystals exerts a very small influence 
on the matrix elements of the exciton-phonon interaction. Note that the exciton-phonon 
coupling due to this mixing in PbS quantum dots calculated using the four-band envelope- 
function formalism M is over 2 orders of magnitude smaller than the value of CdSe quantum 



S 



dots at all radii. Hence, the influence of boundary structure imperfections on the probabil- 
ities of one-phonon optical transitions becomes of paramount importance. This gives us a 
reason to consider electron-hole states in PbS within the model of simple bands. The influ- 
ence of interface imperfections is modeled by the potential expanded in spherical harmonics 

00 1 / r \i 



-1 m=—l 



U imp (r) = J2 E U im y Y lm (0, <p) , (26) 



which obeys the Laplace equation. Exciton quantum states \fj,j) in Eq. (p3[) are calculated 
in the first perturbation approximation on the potential U- imv> (r) . Calculated scattering 
intensities are then averaged on the Gaussian distribution of random amplitudes Z7j m with 
the variance Uq as a fitting parameter. Eq. (BBI) describes, in particular, the electrostatic 
potential, induced by trapped charge. This charge can build up on the nanocrystals during 
the steady-state Raman measurements RJT2]] . 

In Raman scattering experiments |6|-§], the frequency of the incoming light is chosen 
to be in resonance with the exciton ground-state energy. Hence, only several lowest levels 
should be considered when calculating scattering amplitudes (|2~3|) . Two types of selection 
rules manifest themselves in the phonon-assisted Raman scattering. The first group of 
selection rules results from the symmetry properties of the dipole matrix elements d and 
d%: (i) for parallel polarizations (e 7 || e s ), the projection of the exciton angular momentum 
is preserved in the scattering process, M = M'; (ii) for crossed polarizations (e 7 _L e s ), 
M = M' ± 1. 

The second group of selection rules is determined by the symmetry properties of matrix 



elements of the exciton-phonon interaction amplitudes {jij f3 v . Only the following 

phonons take part in the one-phonon scattering: (i) phonons with m = 0, in parallel polar- 
izations; (ii) phonons with m = ±1, in crossed polarizations. The selection rules for the one- 
phonon Raman scattering in spherical nanocrystals were analyzed in Ref. supposing a 
non-degenerate parabolic band structure, and the conclusion drawn was that only s-phonons 
(with I = 0) are active in the dipole scattering processes. This conclusion cannot be applied 
to nanocrystals with a degenerate valence band, where due to band mixing s-phonons and 
(i-phonons (with I = 2) appear to be active in the one-phonon Raman scattering. 

In the multi-phonon scattering amplitudes of Eq. (|23|) with the phonon number K > 2, 
each matrix element of the exciton-phonon interaction contains at least one intermediate 
quantum state, which can be optically active or inactive. Since in the scattering amplitudes 
of Eq. (|23|) , the summation is carried out over all intermediate states, phonons of arbitrary 
angular momenta can participate in multi-phonon transitions assisted by more than one 
phonon. It should be mentioned, that the adiabatic approximation allows participation of 
only s-phonons in multi-phonon scattering processes |J. 

It is worth noting that the one-phonon Raman scattering assisted by p-phonons is for- 
bidden in the case when a nanocrystal has the inversion symmetry. However, the potential 
due to imperfections [Eq. (|26|)] contains both even and odd terms and hence breaks this 
symmetry down. Owing to the breakdown of the inversion symmetry, the p-phonons (and 
other phonons of the odd parity) can be active in the one-phonon Raman scattering. 
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IV. DISCUSSION OF RESULTS 



For the calculation of Raman spectra for CdSe quantum dots, the following values of 
parameters are used: £1(00) = 6.23, £i(0) = 9.56 (Ref. pTCf), the electron band mass m e = 
0.11m , the bulk energy band gap E g = 1.9 eV (Ref. f5"2f), the Luttinger parameters 71 = 
2.04, 72 = 0.58 (Ref. [0), and e 2 = 2.25 (Ref. §). 

For PbS quantum dots, we have taken E g = 0.307 eV (Ref. [0]), E\ (00) = 18.5, 
E\ (0) = 190 (Ref. p5|). The energy of the size quantization in the PbS nanocrystals under 
consideration (with R ~ 1 to 2 nm) is considerably larger than the bulk band gap. For 
such a large energy, the dispersion law of both an electron and a hole is substantially non- 
parabolic, so that we cannot use the effective-mass approach even for simple bands. In this 
view, we determine exciton energies from the experimental absorption spectra of Ref. ||. 
The peak positions observed in Ref. || for PbS nanocrystals of the radius 1.5 nm correspond 
to the energies E exc (ls, Is) « 2.06 eV, E exc (lp, lp) « 3.14 eV, E exc (ld, Id) « 4.18 eV. Since 
the conduction and valence bands in PbS are known to be approximately mirror images 
of each other, the energies of transitions between the lowest states of an electron and of a 
hole can be estimated as E e ^ (lp) — E e (h) (Is) ~ [E e xc (lp, lp) — E exc (ls, Is)} /2 0.54 eV, 
E e{h) (Id) - E e(h) (Is) « [E exc (Id, Id) - E exc (ls, Is)} /2 « 1.06 eV. 

The Raman scattering spectrum, calculated for a Gaussian ensemble of CdSe quantum 
dots in glass and compared with experimental data of Ref. ||, is shown in Fig. 1. The 



relative size dispersion 5 = J (AR 2 )/ (R) is taken 10% (see Ref. ||). The parameter = T 

describing the linewidth of exciton states is taken T = 0.3co> 1)L O) that is close to values 
measured in Ref. |5l| . 



The structure of the phonon-assisted Raman scattering peaks for the spherical CdSe 
nanocrystals appears to be as follows. According to the selection rules discussed above, only 
s- and (i-phonons are active in the one-phonon spectrum. In crossed polarizations, transi- 
tions assisted by s -phonons are forbidden. Hence, the one-phonon spectrum is determined 
only by <i-phonons with m = ±1. In CdSe nanocrystals, the probabilities of transitions are 
determined mainly by the considerable effect of mixing between bands of heavy and light 
holes, while the role of interface imperfections is of minor importance. It is worth noting 
that, as a result of non-adiabatic transitions within the state (lS,lS3/2), the contribution to 
scattering intensity due to cZ-phonons is significantly enhanced compared to that calculated 
|51,5"7fl in the adiabatic approximation. This enhancement is a result of the Jahn- Teller ef- 
fect. In the multi-phonon scattering processes (with K > 2), the participation of phonons of 
other types, in particular, p- and /-phonons, is permitted under the above selection rules. It 
is also worth noting that transitions assisted by p-phonons through excited states (lS,lP3/2) 
and (lS,lP5/2) bring a key contribution into the multi-phonon scattering probabilities for 
spherical CdSe quantum dots with typical radii R ~ 1 to 4 nm. Recently, in Ref. jKJ, the 
multi-phonon Raman scattering in spherical semiconductor nanocrystals has been treated 
within the model |L5]| for the phonon spectrum and for the electron-phonon interaction. In 



Ref. |[46|| , no imperfections and no non-adiabatic effects are taken into account. Relative 
overtone intensities, calculated in Ref. [[16[] using the effective-mass approximation for exci- 
tonic states, quantitatively differ from the experimental data of Ref. . In order to provide 
a quantitative description of the Raman spectra, in Ref. [46] a correction factor has been 
introduced. As seen from Fig. 1, our results compare well with those of Ref. without 
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using any correction factors. 

The Raman scattering spectrum obtained for an ensemble of PbS nanocrystals with a 
relative size dispersion ±4% (Ref. ||) is shown in Fig. 2. This spectrum is calculated with 
Uq = 0.036 eV, which is much less than the typical energy difference between exciton levels. 
Furthermore, the value = T = 15 cm -1 corresponds to the experimentally observed peak 
broadening which is attributed to a finite phonon lifetime (see Ref. |§). The fundamental 
scattering intensity in PbS nanocrystals due to both adiabatic and non-adiabatic transitions 
is determined by trapped charge and is proportional to Uq, while overtone peaks are formed 
mainly as a result of non-adiabatic transitions through the lowest excited states of an exciton. 

The main contribution to both fundamental and overtone bands in PbS quantum dots 
comes from the exciton-phonon interaction with lp- and 2p-phonons. This contrasts both 
with the results of the adiabatic approach, which implies a domination of peaks correspond- 
ing to s-phonons, and with those for CdSe quantum dots, where s- and <i-phonons are 
active in the fundamental scattering. Peak positions given by the adiabatic theory differ 
significantly from those of the experimental data ||. 

As a result of considerable LO phonon dispersion in bulk material, the phonon frequencies 
in PbS quantum dots vary substantially for different modes, i. e. iO\ s = 235 cm -1 , uj\ v = 217 
cm -1 , u 2p = 192 cm -1 , U\d = 229 cm -1 , u>2d = 114 cm -1 . Therefore, the phenomena 
contributing to Raman peaks can be easily identified. In particular, main peaks in both 
fundamental and overtone bands can be confidently assigned to p-phonons. It is evident 
from Fig. 2, that a satisfactory agreement exists between the calculated Raman spectrum 
and the experimental results || with respect to both ratios between spectral intensities and 
peak positions. 
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Fig. 1. Multi-phonon Raman scattering spectrum for the ensemble of CdSe nanocrystals 
with (R) = 2 nm compared with the experimental data of Ref. |J. The thin dashed line 
is the luminescence background. The following parameter values are used: £1(00) = 6.23, 

0.11m , 



£i(0) = 9.56 |nj, m e 
(Ref. §). 



1.9 eV 52 
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2.04, 72 = 0.58 H, and e 2 = 2.25 



Fig. 2. Calculated Raman scattering spectrum (the heavy dashed curve) for the ensemble 
of PbS nanocrystals with (R) = 1.5 nm. The experimental spectrum of Ref. || is shown by 
the solid line. The following parameter values are used: E g = 0.307 eV j53|; S\ (00) = 18.5, 
ei (0) = 190 ~ 
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